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Abstract 

We have examined the role of the BCS pairing mechanism in 
the formation of the magnetic moment and henceforth a spin 
glass (SG) phase by studying a fermionic Sherrington-Kirkpatrick 
model with a local BCS coupling between the fermions. This 
model is obtained by using perturbation theory to trace out the 
conduction electrons degrees of freedom in conventional super- 



conducting alloys. The model is formulated in the path integral 
formalism where the spin operators are represented by bilinear 
combinations of Grassmann fields and it reduces to a single site 
problem that can be solved within the static approximation with 
a replica symmetric Ansatz. We argue that this is a valid pro- 
cedure for values of temperature above the de Almeida-Thouless 
instability line. The phase diagram in the T-g plane, where g is 
the strength of the pairing interaction, for fixed variance J 2 /N of 
the random couplings Jy, exhibits three regions: a normal param- 
agnetic (NP) phase, a spin glass (SG) phase and a pairing (PAIR) 
phase where there is formation of local pairs. The NP and PAIR 
phases are separated by a second order transition line g = g c (T) 
that ends at a tricritical point T 3 = 0.9807 J, g 3 = 5, 8843 J, from 
where it becomes a first order transition line that meets the line 
of second order transitions at T c = 0.9570 J that separates the NP 
and the SG phases. For T <T C the SG phase is separated from 
the PAIR phase by a line of first order transitions. These results 
agree qualitatively with experimental data in GdzThi^Rl^. 



PACS numbers: 05.50. +q, 6460.Cn 
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1. Introduction 



Experimental evidence in cuprate superconductors JT| exhibit a very rich 
phase diagram that includes structural, antiferromagnetic, insulator-metal, 
superconducting and spin glass transitions, that depend strongly on the 
dopant concentration. The coexistence of spin glass ordering and super- 
conductivity has been also observed in conventional superconductors doped 
with magnetic impurities ||. Theoretical studies of conventional spin glass 
superconductors have focused in calculations of the superconducting density 
of states in the presence of localized magnetically ordered impurities. These 
systems are well described by a Hamiltonian where the superconducting 
electrons are represented by a conventional BCS Hamiltonian and they inter- 
act with the localized magnetic impurities via the s-d exchange interaction. 
Theoretical studies of superconductive glass models that describe random 
arrays of Josephson junctions have been performed both for classical |4| and 
recently in a quantum model 0. 

Our motivation in this paper is to study the interplay of the mechanisms 
that lead to spin glass ordering and BCS pair formation in a fermionic Ising 
spin glass model with BCS pairing among localized fermions of opposite spins. 
We argue in the Appendix that this effective Hamiltonian is obtained from 
the model of Ref. || by tracing out the degrees of freedom of the supercon- 
ducting electrons to second order in the s-d exchange interaction, when the 
localized spin operators are represented by bilinear combinations of fermions. 
In this case, besides the known RKKY interaction between localized spins we 
obtain an exchange induced pairing interaction between localized fermions, 
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mediated by the superconducting electrons. This model allow us to investi- 
gate the competitions between frustration and double occupation of the sites 
in a half-filling situation. 

Since the introduction of the Sherrington-Kirkpatrick || (SK) model to 
describe infinite-ranged Ising spin glasses, a vast amount of work was devoted 
to the study of analogous quantum spin glass (QSG) models with different 
and interesting low temperatures properties. 

In an early seminal paper [0 , Bray and Moore used Feynman's functional 
integrals formalism with a fictitious time 0<r</3, j3 = ^, to analyze the 
quantum Heisenberg spin glass model. By using the static approximation to 
evaluate the spin-spin correlation functions, they established the existence 
of a phase transition at finite temperatures. This formalism has been ex- 
tended recently to the study of quantum fluctuations in related spin glass 
models f| ||, [fi|. The authors in Ref. |§ and Ref.|| report on unconven- 
tional time (frequency) behaviour of the correlation functions at T = 0. In 
a remarkable later work [|10j Grempel and Rozenberg found the exact nu- 
merical solution of Bray and Moore's equations |7| for S = 1/2, and they 
demonstrate the existence of an ordered spin glass phase below a finite criti- 
cal temperature . Also the spin-spin correlation function Q(t) is found to be 
roughly constant and equal to its classical value within a range of tempera- 
tures around the critical point, what it seems to justify the use of the static 
Ansatz of Ref. J7J at not very low temperatures in the Heisenberg spin glass. 
Other functional integral techniques have been used earlier to study QSG 
models where the spin operators are represented by bilinear combinations of 
fermionic (anticommuting) Grassmann fields, both in the anisotropic (Ising) 
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[ |TT|j and the isotropic (Heisenberg) limits within a replica symmetric 
(RS) theory. 

The static approximation was used, and it turned out to be exact, in the 
fermionic Ising model, while the fermionic Heisenberg model was solved by 
combining the static approximation for the order parameter with an instan- 
taneous approximation for the retarded susceptibility. 

Recent work [|13[] demonstrated the existence of several characteristic tem- 



peratures in both models, with the de Almeida- Thouless |14| instability oc 



curring at a temperature Ti lower than the spin glass transition Tsg- in the 
isotropic fermionic model [0 there exists still a lower temperature T 2 < T\, 
at which the replica symmetry stability is restored. The region of RS insta- 



bility is characterized by a negative entropy in the anisotropic |Tl| fermionic 
model, while the entropy remains positive in the isotropic model [12] but the 
specific heat changes sign in the RS instability region. 



The anisotropic (Ising) QSG model |LJ], [L5| deserves some special dis 



cussion. In this particular case, the spin operator commutes with the 
particle number operator rij s = or 1, and thus it would not be necessary to 
employ the functional integral formulation since the Hamiltonian is diagonal 
in occupation number operators. However, there still remains an important 
difference between the fermionic and classical SK spin glass: in the quantum 
case the diagonal component of the order parameter in replica space is no 
longer constrained to unity. Consequently, the susceptility x emerges in the 
problem with an important new role and the spin glass order parameter has 



to be determined coupled to x- By adding to the fermionic Ising [H, a 



term that favors BCS pairing, the use of functional integrals becomes nec- 
essary as the Hamiltonian does not commute with the particle occupation 
number operators. 

There is a crucial aspect that characterizes the representation of spin 
operators in Fock space, because there are four quantum states at every site 
, two of them non-magnetic, and the quantum statistics that controls 



number occupation can induce unusual phase transitions. In other words, 
the QSG frustration can be disrupted as long as we have access with equal 
probability to the magnetic and the non-magnetic states at each site. In 
fact, a recent paper |T6j has reported tricritical behaviour in the fermionic 



Ising QSG model, within the static approximation, by varying the electronic 
concentration. This raises the question if the effects that come from the 
relative occupation of magnetic and non-magnetic states can be properly 
exploited, and consequently to produce unusual phase transitions even if 
the average occupation per site is kept constant and equal to one. The 
main difference between Ref. [HI and ours resides in the mechanism that 



controls the magnetic moment formation on the sites. They achieve that by 
varying the electronic concentration while we have a pairing mechanism that 
energetically favors the double occupation and therefore non-magnetic states 
in the half-filling situation. 

This paper is structured as follows: in Sec. 2 we study the model derived 
in the Appendix and find the thermodynamic potential, together with the 
saddle point equations for the order parameters. In Sec. 3 we discuss the 
nature of the phase transitions in the resulting phase diagram in the T-g 
plane, where g is the strength of the pairing interaction, for fixed variance 



of the random couplings Jy. Lowering the temperature, for high values of g, 
there is a line of second order transitions from normal to the pair formation 
phase, that ends at a tricritical point T = T 3 and g = g$, characterized by the 
simultaneous vanishing of the two first coefficients in the Landau expansion 



of the free energy [17j]- From there on, the line becomes one of first order 
transitions until it meets the line of second order spin glass transitions at 
T = T c . For T < T c the first order transition line separates the spin glass 
and pair formation phases. All these results were obtained by using the static 
approximation. As we discussed previously, we expect this to be a justifiable 
Ansatz because all the relevant temperatures are of the order f3J ~ 1 [|7| 
, and our theory is not applicable to very low temperatures due to the de 
Almeida-Thouless instability. 

We reserve Sec. 4 for discussions and comparison with other models and the 
experimental data f2j. 



2. General Formulation 

Conventional spin glass superconductors are well represented by a Hamilto- 
nian where the conduction electrons are described by a BCS Hamiltonian 
and they interact via an effective s-d exchange term with randomly localized 
magnetic impurities ||. We show in the Appendix that, when the localized 
spins are represented in terms of fermions, the degrees of freedom of the su- 
perconducting electrons can be integrated using second order perturbation 
theory in the exchange interaction J s d to give rise to an effective BCS pairing 
interaction among the fermions, besides the very well known RKKY inter- 
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action among the localized spins. In the mean field spirit we are lead to 
consider the following Hamiltonian: 



H = H-nN = -Y i JijSlS* j -iiY i E c U 



ij i s=T4 



4t J2 c h4i c n c n (!) 

i,3 



where the operator £>? is defined as 



s t = c \\ c i\ ~ 4i c n , ( 2 ) 

/i is the chemical potential, c\ s (cj S ) are fermions creation (destruction) op- 
erators and s =| or J, indicates the spin projection. The coupling Jij is an 
independent random variable with the distribution 



P(J y ) = e- J ^Jy N/27rJ 2 . (3) 

The first two terms in the Hamiltonian of equation ([!]) describe a fermionic 
Ising spin glass ]11 , |l^, |15| while the last term is a BCS-like pairing interaction 



and corresponds to the mechanism that favors the double occupation of sites 
0- 

Our ultimate goal is to reduce this problem to a one-site problem. Func- 
tional integration techniques have proved to be a suitable approach for dis- 
ordered quantum-mechanical many-site problems, as it has been for classical 



problems fL9| . Furthermore, this formulation showed to be quite success- 
ful to describe the usual superconductive transition with a BCS coupling 
Pm and in the presence of transition metal impurities [2T|. In that case, the 



particle-hole transformation introduced by the use of Nambu matrices within 
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the static approximation made the problem with a BCS coupling solvable, 
because it becomes a mean field theory in momentum and frequency space . 
As the static approximation is exact for the BCS problem and is also exact 
for the fermionic Ising glass JTl], [L2|] , we expect it to give reliable interpolation 
results here for finite temperatures jUj. Our theory is not valid at very low 
temperatures due to the de Almeida-Thouless instability, then we are not 
concerned with the singular behaviour found at T = in other models M, 0]. 



In the Lagrangian formulation [JTT], [L9], £L[] the partition function is ex- 
pressed as 

Z = J D(0*0)e A (4) 
where the action A is given by 

is aT 

In both expressions </>* s (t) and «/>j S (r) are anticommuting Grassmann vari- 
ables, t is a complex time and /3 the inverse absolute temperature. 

In order to apply the particle-hole transformation within the static ap- 
proximation and to make explicit our central approximation, we work with 
time Fourier transformed quantities. Therefore, the pairing part of the action 
becomes 

Apairing = ^ E E P*(typj(ty (6) 

n ij 

where 

(H(n) = Y,M-w)Mto + w) (7) 
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with Matsubara's frequencies uj = (2m + 1)it and Q = 2miT, (m = 0, ±1, ..). 
In the static approximation, we retain just the term Q = in the sum over 
the frequency Q. Hence we get for A vairing 

AUin 9 = % EEt-H e* ±M 2 ( 8 ) 

p=l iu) 

where we introduced the Nambu matrices 

MM = (<f>Uu) M-")) ^M=\ MUJ \\ (9) 



and the Pauli matrices 



l\ / — i \ /l 

^2 = £3 = ) • (10) 

10/ I % o / V o —l 



The spin part of the action can also be written within the static approx- 
imation as: 



A S g = J2P J iJ S i S j (11) 

where, from eq. (E3), 



£ = (12) 

Finally, the free action is expressed in terms of Nambu matrices 

A ° = E^lM Go 1 ^) t^) ■ (13) 

i 

where the free inverse propagator is 

G 1 (uj) = % u + /i a 3 . (14) 

and the total action can be rebuild as A = A + A s * airing + Asg- We are now 
able to follow the standard procedures to get the configurational averaged 



thermodynamic potential by using the replica formalism 

Sl = -~ lim {Zn)ca - 1 (15) 
(3 n->o n 

where the configurational averaged replicated partition function, after aver- 
aging over Jij, becomes 

Z(n) = (Z n ) ca = J D(P a ,<l> a ) exp[E^» Ctfju) £(oo) 

iaw 



AN 

a p=l ioj 

a/3 i 

The notation ip^{uj) means that a replica index a = l,2,..,n has been as- 
sociated to each matrix element. We introduce replica dependent auxiliary 
fields r) a and g Q/ g to linearize the action in eq. QI6]) , and get 

J r+oo r+oo 



| r+oo r+oo 

where rj a = rj aR - irj a i, K = (j$j2)(j^) and 

Hla P ,Va) = [ ^(0;,0a)eX P [EEf » ^'H g» + 
^ a a; 



2 J* y £Q a 0S?S?] (18) 

a/3 



while the matrix ^ a is defined as 
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We analyze the problem within the replica-symmetric ansatz 

W/3 = 1 Qaa=q + X Va(V* a ) = V(V*) ( 2 0) 

where q is the spin glass order parameter and x is related to the static sus- 



ceptibility by x = «■ The complex order parameter rj gives the number 



of particle-hole pairs of opposite spin at each site, as is obtained extremizing 

drj 



Z(n), that is, solving ^-(Z n ) = in eq. ([17]) and the corresponding equation 



for 7]*. This yields 

v* = H(<M-^)<MM> = ( c n c n) > ( 21 ) 

w 

where the brackets indicate both, a statistical average and average over dis- 
order. 

The sums over a in the spin part of the action produce again quadratic 
terms that can be linearized by introducing new auxiliary fields, with the 
result: 

/+oo p— — f+oo p~T~ 
dz-=[\ dt-=I&z)) n (22) 

I{^z) = J D^^^e^t'^c-'^ti") (23) 

where the matrix G _1 (u;) is given by: 

fiuj + \(z,£)+f3fx Pgr) \ 

Q ! H = 24 
V flgv* tu + X(z,0-^J 

and 

\(z 1 i)=f3jJ2qz + pjJni. (25) 
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In eq. (p3|), the differencial D(4>*,4>) stands for ll^Yls d(J)*(uj)d(J) s ((j) and 
the functional integral over Grassmann variables separates into a product of 
integrals over exponentials of quadratic forms, that can be readily performed 
with the result fl§f : 



ln/(£, z) = £ In | (T 1 |= £ ln[(m; + A(z, £)) 2 - 

^ /?V-/?Vh| 2 ]. (26) 



To perform the frequency sum in eq. (|26|) one should have in mind 
that the Nambu formalism introduces a particle-hole transformation in the 
fermions of spin down. Then from eq. ( jT4|) and eq. (|16"D we have that 

^ = (4 or) - (q 4) = (n T ) + (n;> - 1 , (27) 

and the converging factors in the frequency sums should be adjusted to these 
prescriptions, with the result 

J(f , 2) = cosh(A(^, 0) + cosh(/V) . (28) 

where 



fjl = J,* + g^f . (29) 

Not giving rise to confusion, from now on we write 77 in place of | 77 \. 
Introducing eq. ( |28| ) in eq. ( p2|) and using eq. (fl7|) we finally obtain for the 



thermodynamic potential in eq. (|1^) at the saddle point: 

^ = ^V 2 x(2g + y-)+/W- 
- / Dz ln[e^ J x cosh(pjJ2q z) + cosh(/3/x')] • (30) 
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where Dz = dz e ~I^ . We want, on the average, to insure the half-filling 
situation of one-electron per site, thus fixing fj, — in eqs. (|28D and (|30|). 
The saddle point equations for the order parameters that follow from eq. 
are: 

Y =[ Dz oosHPJV^ z) 

T)? sinh^Jy/^) 

[cosh(P J^z) +e (-P 2j2 V cosh(^)] 2 1 J 

r, = 1 f Dz e{ ~" 2j2 * ] SinH(39V) m) 

1 2 J cosh{(3J^ z) + e(-P 2j2 *) cosh(^) ' 1 J 

The replica symmetric solution described here is unstable at low temper- 



atures, when the de Almeida- Thouless [14| eigenvalue Xat becomes negative. 



The calculation of Xat m this model follows as in a previous work |T3| , with 
the result: 



A 



AT 



2 2 r°o [i + e^^x) cosh(^) cosh^Jy^^)] 2 
f loo [ e H^ 2 x) cosh(/3^) + cosh(/?Jv^ ^)] 4 ' 1 j 



For the entropy we obtain: 

| = ^/5 2 J 2 x(2g + x)-2^ 2 + 

/+oo ~ 9 _ I 

Dz ]n[e p J x cosh(pJJ2q z) + cosh(^)] . (35) 
-oo 

We show in Fig. 5 the behaviour of Xat and S/K as a function of the 
temperature for a value of g > g c . We observe a discontinuity in the derivative 
of the entropy from the normal to pairing phase typical of the second order 
transition. For lower temperatures Xat and S/K become negative due to 
replica symmetry breaking. 
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A detailed discussion of the numerical solutions of the saddle point equa- 
tions, as well as the Landau expansion of the thermodynamic potential in 
eq. ( p0|) in powers of order parameters q and r\ is performed in Sec. 3. 



3. Phase Diagram and Tricritical Point 

The numerical analysis of the equations for the order parameters q, rj and x in- 
equations (0), (p2|) and (|33"D allow us to build a phase diagram (temperature 
versus pairing coupling g) where three regions can be identified (see fig.(l)): 

i) For high T and small g, we get a normal phase with no long range order 
where q = and r] = 0. 

ii) Enhancing the pairing coupling g, one gets a phase transition at g = 
g c {T) where there is a new order corresponding to the spin pairing on the 
sites. In terms of the order parameters, that means f] while q = 0. 

iii) As one lowers the temperature, for g < g c {T), the model exhibits a 
phase transition at T = T c where q starts to grow continously but with r\ 
still equal to zero as shown in fig. (2). The behaviour of the order parameter 
q and the susceptibility % = ft X shows a second order transition from a 
normal phase to a spin glass phase. Actually, that situation has been already 
analysed in [JTTJ] where an expansion of equations fl3ll) and (^) in powers of 



q for r] = gives T c = 0.9570J. 

The nature of the transition line given by the equation g = g c {T) is 
far more complex. If T > T c , the numerical analysis shows that rj grows 
continously from zero as one crosses the transition line (see fig. (3)). This 
result suggests that we get a second order transition. Hovewer, when T < T c , 
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the numerical solution of the order parameter in fig. (4) seems to indicate 
that the transition line becomes first order at some point. To investigate this 



question further we perform a Landau expansion p2| of the thermodynamic 
potential j3Q in eq. ( |30| ) in powers of the two order parameters r\ and q, that 
define the symmetries of the pairing and the spin glass phases, while x is 



taken at the saddle-point value in eq. (^T|). We find it is more convenient to 
start expanding in powers of q and we write from eq. 



k=0 

where xi.1i Vi T) is the solution of the saddle point equation 

E^(W,T) q k = 0. 
k=o U X 



We look for a solution of eq.fl3~7|) also in the form of a series 



(36) 



(37) 



x = Xo + Xi q + X2 q ■ 



(38) 



with the result that Xo * s given by 



where 



and 



^p^fo(v,Xo,T) = Xo -^ = 



D = 1 + e^ 13 " J2 Xo) cosh(Pgr)) 



d d 2 , 

Xl = -[ (^/O (^=2^)~ \v3o? , 

d d 2 1 <9 3 d 2 



(39) 

(40) 

(41) 
(42) 
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Introducing eq. (|38| ) into eq. (^) by expanding the fkS in powers of q, we 
finally obtain after some lengthy calculations the compact result: 

P 2 J 2 xl 



where 



ln ( e /3 J Xo + !) + AlT] 2 + A ^ + A3r] G _ BlQ 2 _ B2g 3 (43) 



B 2 



o 



1 , a , 6 r 1 15 30. 



4 7 4 



/3V 



1 



L 2/? 2 J 2 D 2J 
2 /? 6 7 6 



(44) 



(45) 



Dn = e> 32j2 XoDr. 



eP 2j2 x° + l 



(46) 



First we notice that the correct solution of eq. fl37|) implies the exact 
cancellation of the term linear in q in (3fl. The order parameter rj and q 
minimize and maximize @, respectively, j3fl in eq. ( |43|). We obtain then 
that the normal paramagnetic phase is characterized by A\ > 0, B\ > 0; the 
spin glass phase with q ^ 0, r\ = by A\ > 0, Bi < 0,; and the pairing 
phase with q = 0, r\ ^ by A\ < 0, Bi > 0. Lowering the temperature 
for small values of g, B\ changes sign first at T c = 0.9570J and as B 2 > 
this is a second order transition line from paramagnetic to spin glass phase 



that was analyzed in detail elsewhere [pTTfl . For g > g c (T) and T > T c we 
have Ai < 0, B x > 0, and the line A\ — g — g c {T) = is a second order 
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transition line if A2 > 0. A quick glance at A2 shows that it is positive at high 
temperatures and negative at low temperatures, then we identify the point 
g 3 , T 3 where A\ = A 2 = as a tricritical point, in agreement with the known 
criteria At this point, a line of second order transitions becomes a first 



order transition line. From eq. fl43|) and eq. (|39|) we obtain T 3 = 0.9807 J, 
g 3 = 5.8843J. The expansion of A\ around the tricritical point gives for the 
critical line: 

A l = -0.0566(> - 9-i) + 0.02575(T - T 3 ) = (47) 

and 7] « (Tp-T) 1 / 2 for T 3 < T < T P , where T P = T 3 + 2.2(g-g 3 ). For T <T 3 
and g < g 3 , the transition from the pairing to the spin glass phase becomes a 
first order transition. Tricritical behaviour has been found previously |H| in 



the fermionic Ising spin glass model with charge fluctuation, and a discussion 
of the relation between this model and ours is left for Sec. 4. 

We show in Fig. 5 the behaviour of A at and S/K as a function of the 
temperature for a value of g > g c . We observe a discontinuity in the derivative 
of the entropy from the normal to pairing phase typical of the second order 
transition. For lower temperatures A^r and S/K become negative due to 
replica symmetry breaking. 



4. Conclusions 

We study in this paper the interplay of the mechanisms that leads to spin 
glass ordering and BCS pair formation in a soluble mean field model Hamil- 
tonian for a fermionic quantum spin glass with a BCS pairing between local 
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fermions. As we show in the Appendix, this model would describe the spin 
dynamics of a superconductive spin glass |2|, |3|] and allow us to study the 
role of the pairing mechanism as control of the site occupation and the local 
moment formation. Comparing our results with previous work that exhibit 



tricriticality [15], |16| , we can see that the pairing order parameter rj enters in 
an effective chemical potential // in eq. (^9|), and our equations would reduce 
to theirs if we make g = 0, fi ^ 0. As we are insuring here half- filling on the 
average, we get fi' = gr] and by varying g we change the site occupation by 
favouring doubly occupied states. 

As a result we obtain the phase diagram in fig. (1) where we observe 
a normal paramagnetic (NP) phase at high temperatures with q — r] — 0. 
By lowering the temperature for g < g c {T c ) we encounter a second order 
transition line from the NP phase to the SG (spin glass) phase at T c = 
0.9570J. For g > g c {T c ) the second order transition is from the NP phase 
to the PAIR ( pairing formation) phase if T3 < T < Tp(g). At T = T3, 
g = #3 there is a tricritical point where the pairing transition becomes first 
order. This point is almost indistinguishable in the fig. (1) from the point 
T = T c , g = g c . For T < T c the line g c {T) becomes a first order transition 
line separating the SG and PAIR phases. The phase diagram obtained in fig. 
(1) is in good qualitative agreement with the experimental results of ref. || 
for Gd a; Th 1 _ 3 .RU2 samples, by assuming that the ratio J/g is proportional 
to the Gd concentration. This assumption is reasonable, as we show in 
the Appendix that the effective value of g is proportional to the number of 
superconductive pairs, that for a given temperature decreases drastically with 
the concentration of magnetic impurities, leading to an increase in J/g. To 
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conclude, we studied a model for a fermionic SK spin glass with BCS pairing 
among the local fermions that is soluble by reduction to a one site problem. 
Although this model originates in the description of the spin dynamics of 
conventional spin glass superconductors, we hope that these results may be 
also relevant for the study of strongly correlated fermions systems through the 
localized one site approximations p3| . It is possible to extend the analysis of 
Ref. [1(J for the study of the time correlation functions in the present problem, 
but it will be the subject of a future work. 
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Appendix 



Conventional spin glass superconductors are usually represented by a system 
of conduction electrons with BCS coupling interacting with localized spins 
[H]. Using Gorkov's decoupling scheme the Hamiltonian is: 



Halloy — XX Cfc ~~ ^) a \s a ks ~ Yl^^U^-kl + ^k a -kl a k]} 



ks 



-J sd Si^i (48) 



where Si is the magnetic moment localized at the random site Ri and Sj is 
the local spin density of the conduction electrons 



J(k-k').Ri n \ 

kk' ss' 



where a ss ' indicates the elements of the vector Pauli matrices and a ks (ak s ) 
are the usual creation (annihilation) operators for superconducting electrons. 
The order parameters A\ are to be determined self-consistently from the 
equations of motion, but here we consider them to be phenomenological 
parameters. In this paper we choose to represent the localized moments by 



a bilinear combination of fermion operators [12 



S i = C \s°ss'Cis> (50) 
ss' 

as we did in eq. @. 

The partition function for the superconducting alloy may be written in 
terms of functional integrals as we did in eq. (f|) and eq. (||): 



J alloy 



is ks 
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where the action is now given by 

^ is ^ ks ^~ 

-H all<yy {4>*{r)^{r)}dr (52) 

and we introduced the Grassmann fields (jf ks , 4>k s for the conducting band. 
Using the Nambu formalism as we did in Sec. 2 we may write the partition 
function 



J alloy 



1] Dtfl, fa) e A ° j I] D{<f>t s , <j> ks ) e A »™ +A °* (53) 

is ks 



where A is the action for non-interacting fermions given in eq. (|T^) and 
Abcs is the action for the superconducting electrons: 

Abcs = E^lH ^r 1 ^) tk^) • ( 54 ) 

ku) 

where similarly to eq. © and eq. fl24j) 



/ 0fct(u;) \ 

f fe H = (^ T H 0_ fci (-^)) g» = , J (55) 

V^-fci(-w)/ 

/iw - /3(e fe - fi) (3A k \ 

The s-d exchange part of the action is given by 

A sd = -J^EE&HMW (57) 

i f2 



where from eq. (f49[) and eq. ( |50| ) 

= EE^( w + fi )^Mw) (58) 
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m) = EEE e ' (i "" ) '^L(w + ^)4'^( w ) (59) 

ss' kk' w 

We indicate by uj = (2n + 1)tt and Q = 2mr the fermionic and bosonic 
Matsubara frequencies, respectively. When we are interested in the localized 
spins dynamics, the conduction electrons degrees of freedom in eq. may 
be integrated out to second order perturbation theory in J s ^ to give the result: 



^ alloy 

where 



jY[D((f>i,<p is )e A ° +A eff^^ (60) 



Z Slot/3 ij 

-kpj sd ) 2 Y.Y< w v(msnn)sj(n) + s;(n)s+(-n)] (6i) 

/ n ij 

and the dynamic interactions V^(Q), Wij(Q) are obtained from the cor- 
relation functions {Sf(Q)Sj(—fl)). By performing the frequency sums as 
indicated in Sec. 2 and approximating Afc « A, we obtain for the static part 
with Q = 0: 



Vf(0) = <W E e *<*"^ E n % + ^~ n ^ (62) 
g P fc + 9J - 



= i# E E B W- B ( k + *l (63) 

' V ^ E\k + q)-E\k) 



where 



n(q) = + 

B(k) = J^Ltanhi^^-) (64) 
£(£) = [(e(£)-//) 2 + |A| 2 ]§ 
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The effective interaction in eq. (|6"2]) was calculated to lowest order in 



|A| and it represents the familiar RKKY interaction |24| that is responsible 
for spin glass ordering. The interaction Wij is of a different character and 
it represents the coupling induced by the exchange of superconducting elec- 



trons. The function B(k) is the matrix element of the pairing operator 



cL-cLfci, and since pairing interactions smooth out the jump in the single 



particle occupation number n(k) we can approximate B{k) = 2E ^ at low 



|A| 

temperatures, what gives in eq. (|63D. 

W, «(0) « — n y e^-^ y _ |A|2 =, [-^ ^] (65) 

As the last sum in eq. ( p5|) is very weakly dependent on q it can be approx- 
imated by its value when q = 0, what gives 

1 IAI 2 

Wa (0) « — V ^^<% (66) 
uV ; E\k) 3 1 ^ 

Introducing eq. (|62"D and eq. ( |66"D in eq. ([61] ) we obtain for the static part 
of the interaction 

A? ff (<j>* s ,<f> is ) « -£ J^RKKY)^).^) -sE^(0)Sr(0) (67) 

We argue that the last term in eq. fl6"7|) is the static part of the action 
corresponding to the Hamiltonian H[ = g J2i c lt c «J. c l| c iT that within Gorkov's 
formalism would give rise to terms — g[A\ J2i c i\ c i\ + J2i c !t c 1jJ; what ulti- 
mately justifies our choice of Hamiltonian in eq. ([!]). 
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Figure captions 



Figure 1: Phase diagram as function of temperature and pairing coupling 
g/ J. Solid lines indicate second order transitions while the dotted line indi- 
cates a first order transition. The tricritical point T 3 , g 3 is shown in detail 
in the diagram. The points where \at becomes negative are represented by 
the dashed line. 

Figure 2: Temperature behaviour of q (solid line), and X = f (dashed line) 
for g = 0.5J < g c . Here 77 = 0. 

Figure 3: Dependence of the order parameter rj (dotted line) and the pa- 
rameter x (dashed line) with the coupling g/J for T > T c , where T = 1.5J. 
We show in detail the continuous behaviour of r\ around g c . 
Figure 4: Dependence of the order parameters rj (dotted line), q (solid line), 
and the parameter \ (dashed line) with the coupling g/J for T < T c , where 
T = 0.75 J. At the transition point both 77 and q have discontinous behaviour 
indicating a first order transition. 

Figure 5: The Almeida-Thouless eigenvalue (solid line) and entropy (dashed 
line) for g = 6.2J. 
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